PSEUD OMOMENTS OF THE RIEMANN 
ZETA-FUNCTION AND PSEUDOMAGIC SQUARES 

BRIAN CONREY AND ALEX GAMBURD 

Abstract. We compute integral moments of partial sums of the 
Riemann zeta function on the critical line and obtain an expression 
for the leading coefficient as a product of the standard arithmetic 
factor and a geometric factor. The geometric factor is equal to 
the volume of the convex polytope of substochastic matrices and 
is equal to the leading coefficient in the expression for moments of 
truncated characteristic polynomial of a random unitary matrix. 



1. Introduction 



1.1. Moments of the Riemann zeta function. The Riemann zeta 
function is defined for Re(s) > 1 by 

-i 
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(D «•> = £;? = n 

71=1 p 

As is well-known C nas nieromorphic continuation to the whole 
complex plane with a single simple pole at s — 1 with residue 1. Fur- 
ther, it satisfies a functional equation, relating the value of £(s) and 
the value of £(1 — s), 

(2) C0O=x(*)C(l-a), 
where 

(3) -nss-l. 



Following the standard notation we write s = a + it. 

The problem of computing the moments of £ on the critical line 
a = \ is fundamental, difficult and longstanding. 

The second moment was obtained by Hardy and Littlewood ^H] i n 
1918: 



(4) 



1 

T 



T 



dt ~ logT, 
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the fourth moment was obtained by Ingham [22] in 1926: 



(5) 
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T Jlt 



dt 



2tt : 



log 4 T. 



The asymptotics of higher moments is not known. It has long been 
conjectured that 



(6) 
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T 



2k 



dt ~ c fc log fc T. 



In 1984 Conrey and Ghosh [7] gave the moment conjecture a more pre- 
cise form; namely, they conjectured that there should be a factorization 



gk^k 



(7) Ck= r(i + P) 

where is an arithmetic factor given by 

k 2 oo 



p \ F/ j=0 



pi 



and gu is a geometric factor, which should be an integer. Using Dirichlet 
polynomial techniques Conrey and Ghosh [8] conjectured that = 42 
and Conrey and Gonek |9] conjectured that g\ = 24024. 

1.2. The Riemann zeta-function and characteristic polynomi- 
als of random matrices. In the past few years, following the work 
of Keating and Snaith |23j . Conrey and Farmer Hughes, Keating 
and O'Connell [2011^5 and Conrey, Farmer, Keating, Rubinstein, and 
Snaith [U] it has become clear that the leading order asymptotic of the 
moments of the Riemann zeta function can be conjecturally understood 
in terms of corresponding quantities of the characteristic polynomial of 
the random unitary matrices. Let M be a matrix in U(N) chosen 



itii 



uniformly with respect to Haar measure, denote by e 
eigenvalues, and consider the characteristic polynomial of M: 



D i8 N 



its 



N 



(9) 



P M (z) = det(M-zI) = Y[{ 



e i0 > - z) 



Keating and Snaith (see also Baker and Forrester ^J) computed the 
moments of Pm with respect to Haar measure on U(N) and found that 



(10) 



M N (s) = E { 



<u(N)\Pm{z)\ =11 r(j + s) 2 • 
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They also showed that 

jfrW 0(1 + .)' 

v ; iv-oo 7V S G(l+2s)' 

where G(s) is Barnes double Gamma function satisfying G(l) = 1 and 
G(z + 1) = T(z)G(z). For s = an integer 

(12) G{l + k)2 TT i! 

0(1+24) J-j(j + n)!' 

For = 1,2,3 the quantity above is in agreement with the value of g k in 
the theorems of Hardy and Littlewood, and Ingham and the conjecture 
of Conrey and Ghosh. 

The conjecture of Keating and Snaith (considerably refined and 
extended in is as follows: 

Conjecture 1 (Keating and Snaith |2*3*j). 



(13) \ 



T 



<.|j + « 



2k 



2 



dt ~ a k g k log T 



where a k is an arithmetic factor given by (JHJ) and g k is a "geometric" 
factor (here the notation is different from (J7j)) given by 



(14) g k = Mm ^w\ p »Vn = TT 

1.3. Characteristic polynomials of unitary matrices and magic 
squares. The moments of the secular coefficients of the random uni- 
tary matrices have also been recently investigated. If M is a random 
unitary matrix, following the notation preceding equation @ we write: 

(15) 

N N 

P M (z) = det(M - zl) = J](e^ - z) = (-1)" ^ Sc,(M)^(-iy, 

where Scj(M) is the j-th secular coefficient of the characteristic poly- 
nomial. Note that 

(16) Sci(M) = Tr(M), 
and 

(17) Sc N (M) = det(M). 



\2k k ~ l 
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Moments of the higher secular coefficients were studied by Haake 
and collaborators [13 [TH] who obtained: 

(18) E^Sc^M) = 0, E^ISc^M)! 2 = 1; 

and posed the question of computing the higher moments. The answer 
is given by Theorem ^ which we state below after pausing to give the 
following definition. 

Definition 1. If A is an m by n matrix with nonnegative integer entries 
and with row and column sums 

n 

Ti = CLij, 
3=1 



Cj — 2 J a ij j 
i=l 

then the the row-sum vector row(A) and column-sum vector col (A) are 
defined by 

row(A) = (n, . . . ,r m ), 
col(A) = (ci, . . . ,c n ). 

Given two partitions /i = (/xi, . . . , /i m ) and /2 = (jli, . . . , fi n ) (see |2E] 
for the partition notation) we denote by the number of nonnegative 
integer matrices A with row(A) = fi and col(A) = /L 

For example, for /i = (2, 1, 1) and Ji = (3, 1) we have = 3; and 
the matrices in question are 
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For /i = (2, 2, 1) and \x = (3, 1, 1) we have = 8; and the matrices 
in question are 
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We are ready to state the following Theorem, proved in jTI] 
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Theorem 1 ([O]). (a) Consider a = (a\, . . . , ai) and b = (61, ...,&/) 
with dj , bj nonnegative natural numbers. Then for N > max w2\ jttj , Yl 
we have 

1 

(19) E Un n(S Cj (M))^(Sc,(M)) 6j = 

i=i 

Here fi and jl are partitions \x = (l ai . . . l ai ) , jl = (l bl . . . l bl ) and 

is the number of nonnegative integer matrices A with row(A) = fi and 

col(A) = jl. 

(b) In particular, for N > jk we have 1 

(20) E u(N) \S Cj (M)\ 2k = H k (j), 

where Hk(j) is the number of k x k nonnegative integer matrices with 
each row and column summing up to j - "magic squares". 

1.4. Magic Squares. The reader is likely to have encountered ob- 
jects, which following Ehrhart 14, are referred to as "historical magic 
squares" . These are square matrices of order k, whose entries are non- 
negative integers (1, . . . , k 2 ) and whose rows and columns sum up to 
the same number. The oldest such object, 



(21) 



4 9 2 
3 5 7 
8 1 6 



first appeared in ancient Chinese literature under the name Lo Shu in 
the third millennium BC and repeatedly reappeared in the cabbalistic 
and occult literature in the middle ages. Not knowing ancient Chinese, 
Latin, or Hebrew it is difficult to understand what is "magic" about Lo 
Shu; it is quite easy to understand however why it keeps reappearing: 
there is (modulo reflections) only one historic magic square of order 3. 

Following MacMahon |27j and Stanley what is referred to as 
magic squares in modern combinatorics are square matrices of order 
k, whose entries are nonnegative integers and whose rows and columns 
sum up to the same number j. The number of magic squares of order 
k with row and column sum j, denoted by Hk(j), is of great interest; 
see [TU] and references therein. The first few values are easily obtained: 

(22) H k (l) = k\, 



1 We remark that in ^S] the answer is also obtained in the case N < jk: it is 
related to enumeration of magic squares with certain additional constraints. 
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corresponding to all k by k permutation matrices (this is the k-th 
moment of the traces of powers leading in the work of Diaconis and 
Shahshahani jT5j to the result on the asymptotic normality); 

(23) H x {j) = 1, 

corresponding to 1 x 1 matrix [j] (this is the result of Haake and col- 
laborators given in equation iJTBjl ). We also easily obtain #2(7) — J ' + 1, 
i j -i 



corresponding to 
more involved: 



J -* 



but the value of H 3 (j) is considerably 



(24) ".^('"rn'rHr 

This expression was first obtained by Mac Mahon in 1915 and a simple 
proof was found only a few years ago by M. Bona. The main results 
on H k (j) are given by the following theorems, proved by Stanley and 
Ehrhart (see EEH HE E01 EH ) : 

Theorem (Stanley). H k (j) is a polynomial in j of degree (k — l) 2 , 
having "trivial zeroes" at the negative integers, 

(25) H k (-1) = H k {-2) = ■■■ = H k (-k + 1) = 0, 
and satisfying the following "functional equation": 



(26) H k (-k-j) = (-l) k - 1 H k {j). 

It can be shown that the statements above are equivalent to 

(27) ^H k (j)x> = _ rVfc -i )2+ i > = A: - 3fc + 2, 

i>o 11 XJ 

with /iq + /ii + . . . hd 7^ and hi = hd-i- 
For example, 

ftW) = i/ + |j- + ^- + |i + i. 

and 

1 + x + x 2 



(1 - a;) 5 
j>o v y 



rr , v „■ 1 + 14x + 87x 2 + 148x 3 + 87x 4 + 14a; 5 + x 6 



)x :) 



(l-x) w 
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Theorem (Ehrhart). The leading coefficient of H k (j) is the relative 
volume of B k - the k-th. Birkhoff polytope, i.e. leading coefficient is 

equal to v fy-i k) ■ 

By definition, the k-th Birkhoff polytope is the convex hull of per- 
mutation matrices: 

k k 



(2* 



(xij) g 



Xij > 0; 



x 



8=1 



3=1 



X 



In the example above, volf^) = | x 9. 



1.5. Pseudomoments of the Riemann zeta-function and pseu- 
domagic squares. The purpose of this paper is to prove the following 
result: 

Theorem 2. Let a k be the arithmetic factor given by equation (JHJ). 
Then 

2k 

dt = a klk (\o g xf + O U\ogXf- 1 



(29) 



1 f T X 1 



Here ^ k is the geometric factor, 7^ = vol(V k ), where V k is the convex 
polytope in M k defined by the following inequalities: 



(30) V k 



\X{j) G 



k k ~\ 

^ , ^ X{j ^ 1 , ^ ^ X{j ^ 1 / . 
i=l j=l ) 

The connection with the characteristic polynomials of unitary ma- 
trices is as follows. From Theorem Q it follows that if we consider 
truncated characteristic polynomial 



(31) Pm,i( z ) 

= ^S Cj (M); 

j=0 

we have for N >lk 

(32) ^ u{N) \P M ^z)\ 2k = G k {l), 

where G k {l) denotes the number of k x k nonnegative integer matri- 
ces with row and column sums less than or equal to / (referred to as 
"pseudomagic squares" by Ehrhart [H]) : 
(33) 

k k ~\ 

,j , ^ ^ Xij ^ / , ^ ^ X jj ^ / / . 

1=1 j=l J 



G fc (Z) = card { (ar y ) G Z A 



.r 
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Ehrhart proved that Gk{l) is a polynomial in / of degree k 2 with 
leading coefficient given by jk = vol('Pfc); in fact Gk{l) = card (iVk H Z fc2 
For example, 

G 2 (l) = l(l + l)(l + 2)(l 2 + 3l + 3), 



6 

and 



vol(P 2 ) = -. 

o 



Hence we can rewrite the geometric factor 7^ in a manner similar to 
the expression for in (|T4*|) as follows: 



Eu(Uc)\Pm,i{z, 



\2k 



(34) 7 fc = Km , 

l— >oo < K 

The proof proceeds as follows. In section |2] we obtain an expression 
for 7^ in terms of a multiple complex integral. In section El we express 
the left-hand side of (|29|) as a multiple complex integral and then show 
that the leading terms in the two resulting expressions are equal. 
Acknowledgements. It is a pleasure to thank Dan Bump for pointing 
out that (}3*2"j) is a consequence of Theorem ^ 

2. PSEUDOMAGIC SQUARES 

Let Gk{l) denote the number of k x k nonnegative integer matri- 
ces with row and column sums less than or equal to I given by (J3~3~|) 
(we remark that Hk + i(l), the number of magic squares, is obtained by 
imposing an additional diophantine inequality Yliij x ij — — 1)0 

We have the following expression for Gk(l) as a multiple complex 
integral: 

Proposition 1. Notation being as above we have 
(35) 

Iff (w 1 ...w k z 1 ...z k )~ l - 1 Yl k l=1 dw i U k j= idz j 



( 2 ™p j 7ki=j n . .(1 _ WiZj) n j =l( i _ Wi) n* =l( i _ %: 



The proof follows the approach in |2j, which we now review. 

Let U 1 denote an n-dimensional integer lattice in R" and let V be 
a convex polytope in M n whose vertices are on the lattice Z n (Vk is a 
convex lattice polytope in Z fc2 ). 

Any convex lattice polytope situated in the nonnegative orthant can 
be described as an intersection of finitely many half-spaces: 

(36) ?={x£ M> j Ax < b} , 
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where A is an m x n integer matrix and b e Z m . Consider now the 
function of an integer- valued variable I describing the number of lattice 
points that lie inside the dilated polytope VP: 



(37) 



L(P,l) = card{/P n Z n }; 



with this notation = L(Vk,l)- Denote the columns of A by 

Cj, . . . , c n . Using multivariate generating functions it is proved in [2] 
that for the lattice polytope V given by (j3T?j) we have the following 
expression for L(P, I): 



(38) L(P,l) 



n 



m 

3=1 Z j 



lb~-l 



M= e , nr=i(i-^)nr =1 (i 



■ dz. 



In the expression above we use the standard multivariate notation x y 



Now for Vk the defining system of diophantine inequalities is given 
in ()30|) : the corresponding A is a (2k x k 2 ) matrix given by 



/l 



(39) 



A 



\ 



V 



1 



1/ 



and b = (1,...,1) e Z 2fc . Proposition now follows from ([38)1: for 
notational convenience we have split the variables into two groups 
tui, ...,w k and z x , ...,z k . 



3. Proof of the Theorem 



By the mean-value theorem for Dirichlet polynomials due to Mont- 
gomery and Vaughan [2E] we have 



(40) 



lim — 





x 1 


f 

h 


n=l U ^ +U 



2k 



n=l 



where dk,x{n) is defined by: 

4,x(n) 



h...l k =n 
ll,...,lk<X 
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Consequently we have: 



(41) lim - 





x 1 


f 

/o 


n=l n 



2k 



dt 



Kh<X 



1 

y/h... l k m 1 



l\...l k =m\...m k 



Now we use the discontinuous integral 



1 r +i °° X s , | 0, if < X < 1, 



(42) / — ds . 

27r« J c _ ioo s I 1, if X > 1 



where c > to pick the terms of the Dirichlet series which are less than 
X. Denoting the integral in equation (|42|) by J we can now express 

the right-hand side of (|41|) as follows: 

(43) 

b h 

/f m m m m 

J i= i Ui j=i v i 



(2) (2) 



where 
(44) 

F(u 1 ,...,u k ,v 1 ,...,v k )= ^ 



Ii>l l l "A ™1 



ii-.-ife=mi...mfc 



To simplify notation let u = (ui, . . . , Uk), v = (t>i, . . . , ffc), du = 
dtti . . . duk and c?v = dv\ . . . dvj~- 

Now since for a multiplicative function g(n) we have the Euler prod- 
uct identity: 



oo 

(45) 9( n ) = Hi 1 + 9(P) + ^(P 2 ) + 9{P 3 ) + ■■■), 
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it follows that 
(46) 



F(m, ...,u k ,v 1 ,...,v k ) = Yl 



( 

oo 

E E 

n=l aiH hctfe 

\ a;>0,/%>0 



11 



\ 



1 



„ai(|+«i)H h/3fc(|+^fe) 

n ^ 



p V »j 



+ . . • 



G{u,v)Y[C(^ + u i + v j ), 



where G is an Euler product which is absolutely convergent for \m\ < 



- \v-\ < - 

4' \ u 3\ ^ 4" 

Since 



E 1 = 



QlH hOfc=n 

ai>a,i3j>0 

if we let all Ui and t> ^ be equal to 5 we obtain 

1 oo 

(47) G(*, . . . 8) = JI(1 - -^f E 

and consequently 
(48) 



-2n<5— ra 



n=0 



u limG(u,v)=n(l--) E 



d k (p 



n\2 



pT. 



Summarizing, we have obtained the following expression for the left- 
hand side of 



2k 



( 49 ) T lim hf E^r 

T ^°° T Jo ri2 + 

1 f /"_, vTT^ ySK+^) 

_y...y G ( U) v)nc(i+^+^) T! 



cfu dv. 



Now using the fact that (s — 1)C( S ) is analytic in the entire complex 
plane together with the standard techniques and bounds pertaining to 
C, we obtain that the leading term in (|4"9j) is given by 



(50) 



(2TTi) 2k 

c c 

where we have used (l48l). 



du dv, 
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Write 
(51) 



n w (i 



n 



E 



-Ui—Vj\aij 



A term e _ua e~ v/3 in this expansion is obtained by choosing an N-matrix 
A* = (dijY of finite support with row(A) = a and col (A) = (5. Hence 
the coefficient of e _uct e _v/3 in (|51jl is the number jV^ of N-matrices A 
with row(A) = a and col (A) = (3: 

(52) 

y 1 FT (1 - e~ u ^ v i 

1 ll >3 V a/3 



Further, let I = log X and rewrite the integral 



as follows: 



(53) 



1 

2m 



c+ioo gZs 



0, if I < 0; 

1, if/>0. 



We now express the integral appearing in (JHOj) using (|52~J) and apply 
(J53J) to obtain: 



(54) 



(2m) 2fc 
1 

2fc 



jE(«i+»j) 



(27Tf 



n( da, Jr ' 



du dv 



U; 



^ e -ua e -v/3 



card < (xy) e N A 



a/3 



3=1 



G*(0- 



1=1 



/3<1 

We remark that this proves that the integrals given by (jUj) and 
(|33|) are equal; a direct proof using, for example, a change of variables 
has thus far eluded us. We also remark that the integral expression 
for Gk{l) given by (J35|) has served only as a motivation for the proof 
presented above. We also note that sums related to the expression 
given by the right-hand side of (J4(Jj) were considered in 

4. Generalizations 
Note that in fact we have proved 



rp X ^ 

(55) lim i f dt = a k G k {\ogX) + O ((log X) 



fc 2 -i 
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The proof given in the previous section easily generalizes to yield the 
following result: 

Theorem 3. Let a k be the arithmetic factor given by equation (JSJ). 

Then up to the lower order terms we have: 

(56) 



lim — 



dt ~ a fc C fc (logXi, . . .,\ogX k ). 



Here we assume that Xi = Y m with rrii = 0(1) and Y — > oo; 
(57) 



G k (h,...,l k ) = card< (a; y ) G N A 



X ij — h i X ij — h 

3=1 



8=1 



Finally, we note that in [TT] results analogous to Theorem ^ are 
proved for orthogonal and symplectic group; for example the result for 
symplectic group is as follows: 

Theorem 4. (a) Consider a = (ai, . . . , a;) mtt a, nonnegative natural 
numbers. Let \i be a partition /i = (l ai . . .l ai ). Then for N > Y2ij a j 
and even we have 

i 

(58) E Sp(2N) IJ(S Ci (M)) a ' = NSP^. 

3=1 

Here NSP^ is the number of nonnegative symmetric integer matrices 
A with row(A) = col(A) = \i and with all diagonal entries of A even, 
(b) In particular, for N > jk and jk even we have 

(59) E Sp{2N) S Cj (M) k = S s k p (j), 

where Sl P (j) is the number of k x k symmetric nonnegative integer 
matrices with each row and column summing up to j and all diagonal 
entries even (equivalently, the number of j-regular graphs on k vertices 
with loops and multiple edges). 

We will present analogues of Theorem |2] for L-functions with or- 
thogonal and symplectic symmetries in a forthcoming paper. Here we 
state a representative result for L(s, \d) with Xd{n) — ( ) where d is a 
fundamental discriminant, which has symplectic symmetry. 

Theorem 5. Let b k be the arithmetic factor given by 



, /i , i 
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Then 

ifere is £/te polynomial in I of degree k(k + l)/2 egna/ to £/ie i/ie 

number of k x k symmetric nonnegative integer matrices with row and 
column sums less than or equal to I and all diagonal entries even. 

The connection with the characteristic polynomials of symplectic 
matrices is as follows. From Theorem 0] it follows that if we consider 
truncated characteristic polynomial 

i 

j=0 

we have for N >lk 

from results of Ehrhart it follows that Fk(l) is a polynomial in I of 
degree k(k + l)/2. 
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